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We propose a quantitative criterion to determine whether the coupled quantum systems can 
achieve complete synchronization or phase synchronization in the process of analyzing quantum 
synchronization. Adopting the criterion, we discuss the quantum synchronization effects between 
optomechanical systems and find that the error between the systems and the fluctuation of error 
are sensitive to coupling intensity by calculating the largest Lyapunov exponent of the model and 
quantum fluctuation, respectively. Through taking the appropriate coupling intensity, we can control 
quantum synchronization even under different logical relationship between switches. Finally, we 
simulate the dynamical evolution of the system to verify the quantum synchronization criterion and 
to show the ability of synchronization control. 


I. INTRODUCTION 

In recent years, synchronization effects of two or more interconnected classical systems have aroused comprehensive 
attention because synchronization phenomena are found widely in nature. For examples, Huygens found that two 
clocks with different swings at the initial time will appear synchronization with time evolvement; it was also observed 
that the fireflies glow synchronously and the oscillation of heart cells in human or animal will keep in step with each 
other. At the same time, synchronization effects exhibit inimitable application potential in many fields, such as the 
synchronous transmission of information in the Internet, the synchronous transmission and amplification of signals 
between coupled lasers, the encryption and decryption of signals using chaotic synchronization technology, and so on. 
To date, the synchronization of classical system has gradually become the investigation focus in the numerous fields. 
The groundbreaking work on theoretical exploration of classical synchronization is marked by Yamada and Fujisaka 
who put forward a criterion to judge synchronization behaviors through calculating the Lyapunov exponent of the 
coupling system and obtained synchronization conditions [1]. After this, Pecora and Carroll found synchronization 
phenomena in the electronic circuit and designed a circuit scheme of encrypted communications using synchronization 
techniques, which demonstrates the attractive application prospect of synchronization effects and arouses the intense 
research interest for synchronization theory and application [2]. Recently, many effective synchronization techniques 
have been proposed in order to achieve complete or phase synchronization of classical systems [3-7]. 

Naturally, it is expected to found a similar synchronization phenomenon in quantum systems in order to realize 
the synchronous transmission of quantum information or states due to its unique advantages of synchronization 
effects. However, it is difficult to define precisely some concepts which describe synchronization in quantum systems, 
like “tracks” and “errors”. Even some concepts used in classical dynamics are completely unsuitable to be used 
in quantum dynamics because of remarkable differences of two kinds of systems. Therefore, the relative research of 
quantum synchronization once is thought as unfeasible. The optomechanical system, as a representative of mesoscopic 
systems, has attracted widespread attention and systematic discussion recently [8-11,26]. Mesoscopic systems exhibit 
simultaneously both properties of classical and quantum system under certain conditions because the scale of the 
system is in-between macro-system and micro-system. So some phenomena, no matter what chaos behaviors and 
limit cycle in classical kingdom [12-14,23] or quantum entanglement and quantum coherent in quantum domain [IS¬ 
IS,23], have been observed in optomechanical systems, which provides reliable basis to expand synchronization theory 
from classical to quantum. At 2013, Mari et al. extended the classical synchronization concepts to the quantum system 
[19] and developed quantitative theory of synchronization for continuous variable systems evolving in the quantum 
regime. And in their work, two different measures quantifying the level of synchronization of coupled continuous 
variable are also introduced. Whereafter, some progress has been made in interrelated theories and experiments 
[23-25,27]. 

However in the general case, the synchronization effects are very sensitive to parameters of the systems, such as 
driving field, coupled intensity, and so on. Therefore, it is expected further to investigate and obtain a quantitative 
synchronization criterion in order to determine directly whether the synchronization can be realized. Meanwhile, 
the synchronization criterion can also be regarded as a necessary and sufficient condition of the synchronization 
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effect, which means that the quantum coupling systems can be adjusted and controlled to satisfy the synchronization 
criterion and to realize the aim of quantum synchronization. Further, the controllability and practicability of quantum 
synchronization can be improved. 

In this work, we present a general method for discussing synchronization effects in mesoscopic quantum systems. We 
introduce the first order and second order measurements to describe the expectation value and the fluctuation of error 
respectively and give the necessary conditions to estimate the presence of quantum synchronization effects. Using 
this theory, we design a model based on optomechanical system to realize logic control of quantum synchronization. 
Subsequently, we validate the criterion through the simulation. 

This paper is organized as follows: In Sec.II, the classical synchronization theory is briefly introduced. In Sec.Ill, the 
processing method of quantum mesoscopic synchronization is described and the quantitative criteria for determining 
quantum complete synchronization and quantum phase synchronization are proposed. In Sec.IV, a controllable 
quantum synchronization model base on optomechanical system is designed and the phase synchronization effect is 
discussed. Finally, the summary and the prospects are given in Sec.V. 


II. CLASSICAL SYNCHRONIZATION THEORY 

Considering two classical coupled systems 

dtxi{t) = F{xi{t))+ Ui{xi,X 2 ) 

dtX 2 {t) = F{x 2 {t)) + U 2 {xi,X 2 ) 

where xi{t) and X 2 {t) are state variables of two systems, Ui and U 2 are couplings between systems, respectively. If 
the error X-{t) = |a:i(t) — X2{t)\ —>-0 when t —>■ 00 , the complete synchronization between classical systems is realized. 
If the phases and (j)2{t) of Xi{t) and X2{t) meet ^-(t) = \ t 0, the phase synchronization between 

the systems is obtained. 

The classical synchronization criteria reported previously are mainly to analyze the stability of the error X-(t) or 
phase error (/)_ (t) and to determine whether X- (t) or (/)_ (t) can converge asymptotically to zero through calculating 
the largest Lyapunov exponent. In the next section, we will propose a criterion for quantum synchronization based 
on the classical synchronization theory. 


HI. QUANTUM SYNCHRONIZATION CRITERION 

In the Heisenberg picture, we use quadrature operators qj{t) and Pj(t) to describe two coupled quantum systems 
(here j = 1,2; [qj{t),pj> (t)] = iSjjr). Hence, the error operators q-{t) and p-{t) between the systems can be dehned 
as follows 


q-it) = [qiit)-q2it)]/V2 

p-{i) = [Pi{t)-P 2 {t)\/y /2 

It can be seen from Eq.(2) that the error operators q-{t) and P-{t) are physical quantities which describe the 
differences between the conjugate mechanical quantities of two systems. However, both q-{t) andp_(t) cannot be very 
small simultaneously due to the Heisenberg uncertainty relation. Therefore, it is required to consider synthetically both 
values of two error operators in the quantum synchronization measurement. For this reason, Mari et al. introduced 
the following figure of merit based on Eq.(2) 


S,it) = {q.itf+p.{t)y^ (3) 

Eq.(3) is used to gauge the level of quantum complete synchronization and its value is in the range 0 < 5'c(t) < 1. 

Nevertheless, it is also difficult to discuss quantum synchronization by analyzing Sc(t) directly because it is not 
easy to give a definite criterion for judging the system synchronization or not. On the other hand, the calculation of 
Sc(t) is also complex in quantum systems. Fortunately, the mean value approximation is acceptable in a mesoscopic 
system, meaning that we can write a mesoscopic system operator o{t) in the form of o{t) = {o{t)) +So{t), here (o(t)) is 
the expectation value of the operator at the moment t and it can be regarded as a description of “classical properties”. 
So{t) represents the quantum fluctuation of the operator near its expectation value and the quantum effects of the 
system can be embodied in 6o(t). For a mesoscopic system, So{t) is small but can not be ignored. Then, on the basis 
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of Eq.(2), the error operators of the systems can be rewritten as follow 

= [((ftW) +<5<7 iW) - + 5q2{t))]/= {q_{t)) + 5q-{t) 

P-{i) = [((PiW) - ((P2(0) + 5p2{t))]/y/^ = (p-it)) +Sp-it) 

where Sq-{t) = [dqi{t) — 6 q 2 {t)]/y /2 and 6 p-{t) = [ 6 pi{t) — Sp 2 it)]/y/ 2 .As discussed above, the quantum effects of 
the systems are embodied in 6 o{t). Hence, it is only needed to consider 6 q-{t) and 6 p-{t) in the discussion of the 
synchronization deviation influenced by quantum effects and Eq.(3) can be rewritten as 

S'M = {5q-{tf + 5p.{^X" (5) 

It should be noted that S'^{t) in Eq.(5) is a measurement of the error operator’s quantum fluctuations near its 
expectation value. will equal Sc{t) and it is a synchronization measurement only if {q-(t)) —>■ 0 and {p-(t)) —>■ 0. 

Therefore, we define S'^{t) as a second order measurement to reflect the differences between systems generated by the 
quantum noise even though the classical synchronization conditions are reached. Correspondingly, (g_(t)) —>■ 0 and 
{p-(t)) —>■ 0 can be defined as a hrst order measurement of quantum synchronization. Now that {q-(t)) and {p-(t)) 
satisfy the classical properties, the stability analysis method through calculating the largest Lyapunov exponent of 
the system can be used to determine whether (g-(f)) and (p_(t)) tend stably to zero. Thus, it is only demanded to 
control the stability of {q-{t)) and (p_(f)), that is, we can control mesoscopic quantum synchronization instead of 
analyzing Sc(t) directly. 

Similarly, in the discussion of quantum phase synchronization, the phase error operator can be defined as follow 

(j)- {t) = [<j)i (f) - ( 1 ) 2 (t)]/y /2 (6) 

here (t>j{t) = &i:cian[pj(t)/qj{t)]. We introduce Sp as the measurement of the quantum phase synchronization 

Sp{t) = \{(j).{tf)~" (7) 

Similar to the discussion of complete quantum synchronization, let (j)-{t) = {(j)-(t)) + 5(j)-{t) and rewrite Eq.(7): 

S'p{t) = \{5(j).{tf) ' (8) 

where Sp{t) should be a second order measurement as well under {(f>-{t)) —>■ 0. 

In summary, the synchronization effects in mesoscopic quantum systems can be discussed through the following 
steps: 

a. Write the operator equations of system’s conjugate mechanical quantities in the Heisenberg picture, define the 
error operators and take them as the form of fluctuations near their expectation value, that is, o{t) = {o{t)) + So{t). 

b. Make stability analysis for (o(t)) and calculate the largest Lyapunov exponent of the error equations. If the largest 
Lyapunov exponent is less than zero, the evolution of {o{t)) can tend to zero stably after a certain time, whereas it 
may be ruleless oscillation. 

c. If the largest Lyapunov exponent is less than zero, the following work is to discuss the magnitude of the noise 
iSo{t)) and to calculate 5''(t) and Sp{t) base on Eq.(5) and Eq.(8), respectively. Oppositely, if and Sp{t) keeps 
a constant but not zero, the synchronization between the quantum systems is achieved. 


IV. DESIGN OF THE CONTROLLED QUANTUM SYNCHRONIZATION MODEL AND QUANTUM 

PHASE SYNCHRONIZATION 

A controlled quantum synchronization model is designed based on the quantum optomechanical system in order to 
check the validity of the above-mentioned quantitative criteria. In this model, we can realize quantum synchronization 
control through different logical relationship of the switches, shown in Fig.I. 

Two coupled optomechanical systems are driven by laser and interact mutually through a phonon tunnel and a 
fiber which can be controlled by the open or close of the switches Ki and K 2 . The Hamiltonian of the system can be 
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FIG. 1: Diagrammatic sketch of controlled synchronization model 


given directly after a rotating approximation [ 20 , 21 ] {h = 1 ). 


H = - ga]aj{h] + bj) + iE{a] - Qj)] 

i=i,2 

— n{bibl + b\b2) + A(aia2 + a|a2) 


(9) 


here and Uj are the optical creation and annihilation operators for the system j, and bj are the mechanical creation 
and annihilation operators. Aj and LUj are the optical detunings and the mechanical frequencies, respectively, g is 
the optomechanical coupling constant and E is the laser intensity which drives the optical cavities, /i is the intensity 
of the phonon tunnel and A is coupling constant of the fiber.The switches Ki and K 2 can change /i and A values from 
zero to a positive constant by on and off. After considering the dissipative effects, the following quantum Langevin 
equations can be written in Heisenberg picture through the input-output properties 


dtcii = [—K + iAi + ig(b\ -b 6 i)]ai + E — iXa 2 + V^i 


Ka, 


dtO ,2 = [—K + iA 2 + ig{b\ + b 2 )]a 2 + E — iXai + y/^a 
dtbi = [—7 — iuji]bi + iga\ai + ig,b 2 + y^b™ 

= \~1 ~ + iga ^2^'2 + + \/ 2762” 


( 10 ) 


here k and 7 are the optical and mechanical damping rates, a*" and 6 ™ are the input bath operators, which satisfy 
-b = 5jy5{t - t') and -b b^{t')b^{t)'’) = { 2 nb + l)5jy5{t - t'), where 

nb = - 1 ]“^ 

Using the expectation value and quantum fluctuation to replace operators from Eq.(lO), we can get following two 
equations. 

The parts of expectation value are 


dtA\ — [—K “b iAi -b ig{B* -b E — iXA 2 

dtA 2 = [—K “b iA 2 “b ig{B 2 -b B 2 )]A 2 E — iXAi 
dtBi = [—7 — iuJi]Bi + igAlAi + ifj,B 2 
dtB 2 = [—7 — iuj 2 ]B 2 -b 13 ^ 2^2 -b ig,Bi 


( 11 ) 


where Aj = (uj) and Bj = (bj). 

The parts of quantum fluctuation are 


dtSai = [—K + iAi + ig{Bl + Bi)]Sai + igAi(Sbl + 6bi) — iXSa2 + \P^a 
dtba 2 = [—K + f A 2 -b ig{B 2 + B 2 )]Sa 2 + igA 2 { 6 b\ + ^ 62 ) — iX5ai + y/^a. 
dtSbi = [—7 — iuii]Sbi + igA^Sai -b igAiSa\ + ig,Sb 2 + y/^bi 
dtSb2 = [—7 — iuj2]Sb2 + igA25a2 + igA26a\ + + y/^b- 




bT 


(12) 
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The dynamic properties of the cavity and oscillator can be described by their own conjugate mechanical quantities, 


i.e. 

Xj = (a] + aj)/V2 

Vj = *(«] “ 

Qj = + bj)/V 2 

Pj = Kb] - bj)/V 2 

Substituting Eq.(13) into Eq.(12), Eq.(12) can be expressed in the matrix form 

dtu = Su + ^ 

where m is a vector { 6 x 1 , 6 yi, 6 x 2 , Sy 2 , Sqi, 6 pi, Sq 2 , 6 p 2 )^ and ^ means a input vector 
{ 6 x™, 6 y™, 6 x^^, 6 y^, 6 q'^^, 6 p^{^, 6 q^, 6 p^^)^. S' is 8 x 8 time-dependent matrix. 

! —K — Ai — 2gRe[i3i] 0 A —2gIm[Ri] 

Ai -I- 2oRe[i3il —k —A 0 2gRe[Ril 

0 A -K -A 2 - 2gRe[R2] 0 

—A 0 A 2 -I- 2gRe[R2] —k 0 

0 0 0 0-7 

2gRe[2li] 2gIm[Ali] 0 0 —wi 

0 0 0 0 0 

Vo 0 2gRe[A2] 2gIm[A2] y 

For the convenience of calculation, a covariance matrix C is defined as 

Cij{t) = Cji{t) = -{Ut{t)Uj{t) +Uj{t)Ui{t)) 

and the evolution of matrix C can be determined by Eq.(17): 

dtC = SC + CS^ + N 

where iV is a diagonal noise correlation matrix defined by: 

N^JS{t -1') = +ut')ut)) 


(13) 


(14) 


0 0 0 \ 
0 0 0 
0 - 2 glm[yl 2 ] 0 
0 2gRe[A2] 0 

uji 0 —y 

-7 g 0 

-y —7 UJ 2 
0 —UJ 2 — 7 / 


(15) 


(16) 


(17) 


(18) 


At this point, the dynamical analysis of the model we used has already been finished and the system will evolution 
according to Eq.(ll) and Eq.(14). Subsequently, we are going to discuss the synchronization effects between the 
systems by using those equations. 

We will discuss the phase synchronization between the oscillators of the systems. Firstly, define the “classical” part 
(expectation value) of phase error on the basis of Eq. ( 6 ) [22]. 


where 


(</)_(t)) = liMt)) - {Mm/V 2 = 0 {t)/V 2 


= arg[Rl(t)] = arctan(lm[Rl(t)]/Re[i?l(t)]) 
{4>2{t)) = arg[R2(t)] = arctan(lm[il2(t)]/Re[i?2(t)]) 


(19) 


( 20 ) 


The evolution of 9{t) can be simulated numerically with simultaneous equations (11), (19) and (20) under the 
certain initial conditions and the largest Lyapunov exponent of the error can be calculated by following equation 


L„ = lim — In 

t—>-OC t 


69{t) 


S9{0) 


( 21 ) 


(50(0) and S9{t) in Eq.21 mean the disturbances of the phase errors when t = 0 and t = t, respectively, and Ly can 
be seen as the eigenvalue of Jacobian matrix corresponding to 69. Because the phase synchronization is controlled by 
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the intensity of the phonon tunnel and the coupling constant of the fiber, we calculate the largest Lyapunov exponent 
of phase error with the variation of ^ and A. In calculation, the damping rates and the intensity of driving field 
are assumed to be equal in both systems, but there are the differences in frequencies and initial conditions. The 
values of LO,g,K and 7 are taken as the same as Mari’s work so that the conclusion is more easy to be verified by the 
experiment. Moreover, we properly reduce the intensity of the driving field in order to highlight the coupling function 
in synchronization. Otherwise, too strong driving field will dilute the coupling effect, which leads the system in the 
“forced” synchronous effect. It can be known from Eq.9 and Eq.II that the phonon coupling can directly influence 
mechanical oscillators, however, the photon coupling can only influence them indirectly by changing the light field 
in the cavities. Eor making two switches have similar ability to control synchronization, we reduce the intensity of 
phonon channel(/i) and increase coupling constant of the fiber. Therefore, we calculate the Lyapunov exponent in the 
region of A S [0, 0.2], g G [0,0.01] and the calculation result is shown in Eig.2. 



FIG. 2: Largest Lyapunov exponent of phase error 6{t) with /r and A. Here tui = 1 is the nnit of frequency, and LJ 2 ~ 1.005, 
g — 0.004, K = 0.15, 7 = 0.005, = ujj, rib = 0. 

Instead of analyzing the largest Lyapunov exponent concretely, as it is discussed above, we can determine whether 
the system is in the synchronization state only by comparing the largest Lyapunov exponent with zero. Therefore, 
Fig.2 is redrawn in this form: some parameter regions where the largest Lyapunov exponent is greater than zero are 
projected and marked in red, contrary, the regions where the largest Lyapunov exponent is less than zero are marked 
in blue. Moreover, we also plot a curve of the Lyapunov exponent with a fixed A in order to display it more clearly, 
shown in Fig.3. 



FIG. 3: (a): Gomparison between the largest Lyapunov exponent and zero. Blue areas mean the largest Lyapunov exponent is 

less than zero and red areas mean the largest Lyapunov exponent is greater than zero, (b): Evolution of the largest Lyapunov 
exponent with varied g. Here we set A = 0.16 and other parameters are same with Fig.2. 

If g and A are in the blue regions of Fig.3(a), the largest Lyapunov exponent is less than zero, indicating the 
evolution of the phase error tends to zero stably after a certain time and the systems reach the synchronization. By 
contrast, the largest Lyapunov exponent is greater than zero while /r and A are in the red regions and the systems are 
not synchronous because the phase error tends to random oscillations. Therefore, it can be seen from Fig.3(a) that 
the systems will not synchronize when two systems are not coupled(/r = 0 and A=0). Once there is coupling between 
systems, however, the red regions will be replaced gradually by blue area with the increasing of g and A. We expect 
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to control the synchronization by the switches Ki and K 2 together, i.e. it will be happened only if two switches 
meet the logic “AND”. Then we choose the parameters according to the following principles: the largest Lyapunov 
exponent is greater than zero when /i = 0 and A 7 ^ 0 as well as /r 7 ^ 0 and A = 0, but it must be less than zero as 
/r 7 ^ 0 and A 7 ^ 0 at the same time. According to this, the appropriate parameters can be found in A S [0.14,0.2] and 
/r S [0.004,0.007] recurring to the help of Fig.3 and within the range the systems will achieve phase synchronization 
only if switches Ki and K 2 are closed synchronously. (It is worth to note that alone point with different color from 
ones around can be ignored as an error.) In other words, we can control synchronization effect with switches Ki and 
K 2 by using above characteristics. 

In the above discussion, we give a range of parameters /i and A but not exact values. This is because the “clas¬ 
sical” error can tend to zero for all parameters in that range. Whereas for quantum synchronization, not only the 
expectation value of errors tending to zero is necessary, but also the error fluctuations as small as possible. For the 
sake of reaching a perfect synchronization, we need to select appropriate values of ^ and A and to ensure a minimum 
quantum fluctuations. It is the reason why we calculate the second order measurement Sp{t) of the quantum phase 
synchronization. 

For the calculating of Sp{t) , the matrix C defined in Eq.(16) is transformed firstly as 

C'{t) = C/(t)C(t)C/(t)1' (22) 

here U{t) = ...] and (f)ai = arg(ai(t)), (j)a 2 = arg(a 2 (t)),... . C'[t) can be obtained by substituting 

the matrix C into Eq.(17) and Sp{t) can be expressed as 

S'p{t) = + 6p'l - 2 V 1 V 2 )) ' 

2 2 2 (23) 

Time-averaged Sp{t) is further calculated in order to show directly the size of quantum fluctuation under the 
different parameters. 


S' = lim — 

P T->oo T 


S'{t)dt 


(24) 


and the calculation result is shown in Fig.4. 



FIG. 4: S'p with varied p. and A. In this calculation, we let T = 2000 and other parameters are same with Fig.2 

Although the phase expectation values of every system tend to be equal, the quantum fluctuations between the 
systems under different parameters still influence the perfection of quantum phase synchronization, as shown in Fig.4. 
The fluctuation of system error will be reduced to minimum extent while taking p = 0.004 and A = 0.16, which draws 
the conclusion that the best effect of synchronization has been reached. 











The dynamical evolution of the system is simulated here to test the validity of our criterion. Before the simulation 

TT 

we let /r = 0.004 and A = 0.16 and assume that the initial phase error between the systems is 0(0) = —. The remaining 

parameters are same as ones used in Fig.2. The simulation results are illustrated in Fig.5 in which the unit of ordinate 
is TT. 



FIG. 5: Evolution of the phase error between the systems, (a): Ki and K 2 are both opened = 0,A = 0); (b): Ki is opened 
and K 2 closed (p = 0,A = 0.16); (c): Ki is closed and K 2 opened (/i = 0.004,A = 0) and (d): Ki and K 2 are both closed 
ifi = 0.004,A = 0.16). 


It can be known from Figs.5(a)~(c) that the synchronization between systems will not be achieved as long as any 
a switch of Ki and K 2 is opened. Upon further inspection, we notice that two different systems will never achieve 
phase synchronization in the other parameters {uj,g,K and 7 ) but the same E when couples disappear. Only when 
two switches are both closed, the systems can realize synchronization, shown as Fig.5(d). This result is identical 
with our analysis and it also verifies the quantum synchronous criterion proposed in our work. By the way, the 
logical relationship “AND” of two switches is taken as an example in Fig.3, however, the logical relationships “OR” 
or “exclusive-OR” between the switches can also be selected to realize the quantum synchronization of the systems 
by adjusting appropriate parameters of p, and A . 


V. CONCLUSIONS 

In this paper, we investigate quantum synchronization effects and present the quantitative criteria of complete 
synchronization and phase synchronization between quantum systems. Further, we realize the quantum phase syn¬ 
chronization between coupled optomechanical systems by using our criterisa. Through calculating the largest Lya¬ 
punov exponent, we find that the systems will not reach synchronization unless switches Ki and K 2 are closed 
synchronously (satisfying the logic relation “AND”) when the parameter values are taken at the ranges A G [0.14, 0.2] 
and fjL S [0.004, 0.007]. At the same time, we obtain the information that the fluctuation of system error will reduce 
to minimum while /i = 0.004 and A = 0.16 by calculating the second order measurement Sp{t) of the quantum phase 
synchronization. Finally, the dynamical evolution of the system is simulated in order to test the validity of our 
criterion under above parameters. Since the concrete quantum synchronization criteria have been proposed and the 
control theory of quantum synchronization effects is simple and efficient in the work, other designers can set different 
synchronization conditions to satisfy themselves aims. We believe that our work can bring certain application values 
in quantum communication, quantum control and quantum logical gates. 


Acknowledgements 

This research was supported by the National Natural Science Foundation of China (Grant No 11175033) and the 
Fundamental Research Funds for the Central Universities (DUT13LK05). 















9 


References 


[llYamada T, Fuiisaka H. Stability theory of synchronized motion in coupled- oscillator systems 11. Prog. Theor. 
Phys., 1983, 70(5):1240-1248. 

[2] Pecora L M, Carroll T L. Synchronization in chaotic systems. Phys. Rev. Lett., 1990, 64(8):821-824. 

[3] Selivanov A A, Lehnert J, Dahms T, Hvel P, Fradkov A L, Schll E. Adaptive synchronization in delay-coupled 
networks of Stuart-Landau oscillators. Phys. Rev. E, 2012, 85(l):016201-8. 

[4] Frthauer S, Ramaswamy S. Phase-synchronized state of oriented active fluids. Phys. Rev. Lett., 2013, 
lll(23):238102-5. 

[5] Paredes G, Alvarez-Llamoza O, Cosenza M G. Global interactions, information flow, and chaos synchronization. 
Phys. Rev. E, 2013, 88(4):042920-8. 

[6] Um J, Hong H, Park H. Nature of synchronization transitions in random networks of coupled oscillators. Phys. 
Rev. E, 2014, 89(l):012810-8. 

[7] Kawamura Y, Nakao H. Noise-induced synchronization of oscillatory convection and its optimization. Phys. Rev. 
E, 2014, 89(1):012912-13. 

[8] Kippenberg T J, Vahala K J. Cavity optomechanics: back-action at the mesoscale. Science, 2008, 321(5893): 
1172-1176. 

[9] Thompson J D, Zwickl B M, Jayich A M, Marquardt F, Girvin S M,Harris J G E. Strong dispersive conpling of a 
high-finesse cavity to a micromechanical membrane. Nature, 2008, 452(7183): 72-75. 

[lOjBrennecke F, Donner T, Ritter S, Bourdel T, Khl M, Esslinger T. Gavity QED with a BoseEinstein condensate. 
Nature, 2007, 450(7167): 268-271. 

[lljLin Q, Rosenberg J, Jiang X, K J Vahala, O Painter. Mechanical oscillation and cooling actuated by the optical 
gradient force. Phys. Rev. Lett., 2009, 103(10): 103601-4. 

[12] Zhang K, Ghen W, Bhattacharya M, Meystre P. Hamiltonian chaos in a coupled BEGoptomechanical-cavity 
system. Phys. Rev. A, 2010, 81(1): 013802-6. 

[13] Larson J, Horsdal M. Photonic Josephson effect, phase transitions, and chaos in optomechanical systems. Phys. 
Rev. A, 2011, 84(2): 021804-4. 

[14] Heinrich G, Harris J G E, Marquardt F. Photon shuttle: Landau-Zener-Stckelberg dynamics in an optomechanical 
system. Phys. Rev. A, 2010, 81(1): 011801-4. 

[15] Vitali D, Gigan S, Ferreira A, Bhm H R, Tombesi P, Guerreiro A, Vedral V, Zeilinger A, Aspelmeyer M. Optome¬ 
chanical entanglement between a movable mirror and a cavity field. Phys. Rev. Lett., 2007, 98(3): 030405-4. 

[16] Miao H, Danilishin S, Ghen Y. Universal quantum entanglement between an oscillator and continuous fields. Phys. 
Rev. A, 2010, 81(5): 052307-4. 

[17lHofer S G, Vasilyev D V, Aspelmeyer M, Hammerer K. Time-Gontinuous Bell Measurements. Phys. Rev. Lett., 
2013, 111(17): 170404-6. 

[18] Wang Y D, Glerk A A. Reservoir-engineered entanglement in optomechanical systems. Phys. Rev .Lett., 2013, 
110(25): 253601-5. 

[19] Mari A, Farace A, Didier N, Giovannetti V, Fazio R. Measures of quantum synchronization in continuous variable 
systems. Phys. Rev. Lett., 2013, 111(10): 103605-5. 

[20] Farace A, Giovannetti V. Enhancing quantum effects via periodic modulations in optomechanical systems. Phys. 
Rev. A, 2012, 86(1): 013820-12. 

[21] Mari A, Eisert J. Gently modulating optomechanical systems. Phys. Rev. Lett., 2009, 103(21): 213603-4. 

[22] ln Eq.(2) and Eq.(6), we make the error operators divided by artificially in order to ensure 0 < Sc{t) < 

1 .Discussing the expectation value of the error operator, however, the physical significance of error will not clear 
if make the error operators divided by v^. So, we substitute 9{t) = ((/i'i(t)) — into Eq.(21), instead of 

{cl^.{t)) = {{Mt))-{Mt)))/V 2 . 

[23] Lee T E ,Ghan G-K ,Wang S S. Entanglement tongue and quantum synchronization of disordered oscillators. 
Phys. Rev. E, 2014, 89(2): 022913-10 

[24] Matheny M H, Grau M, Villanueva L G, Karabalin R B, Gross M G, Roukes M L. Phase Synchronization of Two 
Anharmonic Nanomechanical Oscillators. Phys. Rev. Lett., 2014, 112(1): 014101-5 

[25] Manipatruni S,Welderhecker G,Lipson M. Long-range synchronization of optomechanical structures. Quantum 
Electronics and Laser Science Gonference. Optical Society of America, 2011: QWll 

[26] Aspelmeyer M, Kippenberg T J, Marquardt F. Gavity optomechanics, arXiv:1303.0733vl 

[27] Heinrich G,Ludwig M, Qian J, Kubala B, Marquardt F. Phys. Rev. Lett., 2011,107(4):043603-4 


